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We investigate quantised maps of the torus whose classical analogues are ergodic but not 
mixing. Their quantum spectral statistics shows non-generic behaviour, i. e. it does not fol- 
low random matrix theory (RMT). By coupling the map to a spin 1/2, which corresponds to 
changing the quantisation without altering the classical limit of the dynamics on the torus, 
we numerically observe a transition to RMT statistics. The results are interpreted in terms of 
semiclassical trace formulae for the maps with and without spin respectively. We thus have 
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spectral statistics. We also discuss the analogous situation for an almost integrable map, where 
we compare to Semi-Poissonian statistics. 
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1. Introduction 



A long standing problem in quantum chaology is the precise formulation of the conjecture of 
Bohigas, Giannoni and Schmit (BGS) Q which states that, generically, the quantum spectral 
statistics of systems whose classical analogue is chaotic can be described by the average statis- 
tical behaviour of the eigenvalues of large random matrices. In contrast, Berry and Tabor Q 
conjectured that for integrable systems the statistical distribution of eigenvalues is that of a 
Poisson process. 

Whereas in the latter case integrability is precisely defined by the Liouville- Arnold theo- 
rem there is no general consensus as to what should be sufficient conditions for the BGS 
conjecture to hold. It is generally believed that hyperbolicity, i. e. positive Lyapunov expo- 
nents, is sufficient for observing statistics as in random matrix theory (RMT). However, since 
this is a very strong property, one might ask whether weaker conditions would suffice. Actu- 
ally, this question has already been raised by Bohigas, Giannoni and Schmit who originally 
formulated their conjecture for K-systems. The weakest possible property in the hierarchy of 
chaos is ergodicity which, roughly speaking, means that almost everywhere in phase space time 
averages are equal to a phase space average. It is sometimes argued that ergodicity is too weak 
for implying RMT statistics but that instead one should also require the decay of correlations, 
i. e. mixing or even the K-property. 

In order to shed some light on this problem we investigate parabolic maps of the two- 
torus, namely Kronecker and skew maps, which are (uniquely) ergodic but satisfy no stronger 
condition of chaoticity. In particular, they are not mixing, i. e. in general correlations do not 
decay, and have zero topological entropy and hence zero Kolmogorov-Sinai entropy. Thus, 
these maps serve as good examples for soft chaos. 

The quantised versions of these maps are known to show non-generic spectral statistics. 
Even worse, it has been proven Q that some statistical functions do not exist in the semiclas- 
sical limit. Although this might be considered as evidence for ergodicity not being a sufficient 
condition for RMT statistics, we argue that the non-generic behaviour is due to degeneracies in 
both the quantum spectrum and the spectrum of classical periodic orbits. These degeneracies 
can be lifted without changing the characteristic properties of the classical system by perturb- 
ing the quantum system in higher semiclassical order, namely by coupling the quantum map 
to a spin 1/2. We show that in this way one can construct a family of quantum maps which all 
have a non- mixing classical limit but numerically show RMT statistics. We analyse the lifting 
of degeneracies both numerically and in terms of semiclassical trace formulae. Thus we provide 
evidence that generically (which should at least imply the absence of systematic degeneracies 
both quantum mechanically and semiclassically) ergodicity might be a strong enough condition 
for implying the BGS conjecture. 

For skew maps we also investigate an almost integrable situation. The level spacing dis- 
tribution for the almost integrable skew map with spin is close to Semi-Poissonian behaviour. 
This numerical finding is supported by a semi-heuristic analytical calculation of the value of 
the form factor K{t) at r = 0. 

The paper is organised as follows. First we introduce classical and quantised Kronecker 
maps (section § and ^ respectively), which are irrational translations on the two-torus and 
thus represent the simplest possible ergodic system to be studied classically and quantum 
mechanically. In section ^ we discuss their non-generic spectral statistics and relate it to 
degeneracies of classical periodic orbits using a trace formula. The following sections are 
devoted to the study of Kronecker maps with spin, which show generic spectral statistics 



2 



(section ||). We also discuss how a higher order quantum perturbation lifts the degeneracy 
of periodic orbits in a semiclassical trace formula (section ^) . Then in section |^ we turn our 
attention to the slightly more complicated skew maps. Quantum skew maps also show non- 
generic spectral statistics. But again we observe a change to RMT statistics if we add a spin 
contribution. Almost integrable skew maps are discussed in section |8|. We conclude with some 
remarks in section ^. The proof of the Egorov property for the Kronecker map and the trace 
formula for the skew map are discussed in two appendices. 



2. Kronecker maps 

The simplest example for parabolic maps are translations on the two-torus given by 

•p 



K^B : T2 — > t2 



P + P\ , ... (2.1) 

with a, (3 £ S^. If 1, a and /? are linearly independent over Z we have irrational translations on 

the two-torus which are also called Kronecker maps. In these cases we get uniquely ergodic, 

but non-mixing maps (see, e.g., ||5[^]). For a,/3 S Q the maps are periodic. In all other cases 
the torus splits into families of invariant curves. 

For later reference it will be important to study the symmetries of Kap. There are several 
symplectic symmetries vr fulfilling 

vr o Kaj3 o TT^^ = Kaj3 ■ (2.2) 

First of all we have invariance under inversion 



(mod 1) . (2.3) 



Furthermore, there are two continuous families of symmetry operations: the maps commute 
with both translations in position 



and translations in momentum 

W):Q^(^ + '^) (modi), 6eS'. (2.5) 

On the other hand, there is no anti-symplectic symmetry r with 

roKo^f^o = , (2.6) 

and thus Kronecker maps are not invariant under time reversal. 

When discussing the quantisation of the map in the next section, we will also need a 
generating function of the classical map. Due to the independence of the dynamics in position 
and momentum, see (p.l[), one can only give a generating function in a mixed representation. 



Gal3 iqn+l,Pn) = PnQn+l + PQu+I " aPn , (2.7) 

with 

dGa(3{qn+l,Pn) , dGa/S {qn+l,Pn) „s 

Pn+1 = ^ and qn = ^ . (2.8) 

Oqn+l OPn 
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3. Quantum mechanics of Kronecker maps 

Let us first review the basics of quantum mechanics on the two-torus T^, see, e.g., 
Phase space translation operators T{q,p) yield a representation of the Heisenberg group, 

Tiq,p)T{q',p')=e^^'^'P-'^P'^Tiq + q',p + p') . (3.1) 

Because of the topology of the torus, quantum states have to be periodic up to a phase in 
both position and momentum, i. e. 

T(1,0)|V') =e-2^''^i|^) and T(0, 1)1-0) = e^^^''^ , (3_2) 

For Kronecker maps arbitrary ki and K2 are allowed, but the spectrum is not changed by 
choosing different values. Hence for simplicity we choose ki = K2 = 0. The conditions ( ^.2| ) 
require h = {2'n'N)~^ with G IN. Moreover one concludes that a quantum state in position 
(or in momentum) representation is given by a superposition of ^-distributions 

with Ian '■= I1/NI1. Position states are defined by 

{q\Q):=^^6(q-^-m) , with Q G Zat , (3.4) 



and momentum states are given by discrete Fourier transformation of the position states with 
(Q|P) = ^e^^«. (3.5) 



We remark that now the action of translations in either momentum or position is given by 
T(^,0)|Q) = |Q + n) , T(0,f)|Q)=e'^™Q|g) and 



T(f,0)|P)=e-^"^|P) , T(0,^)|P) = |P + 



(3.6) 

m) , 



n,m ^ 1i]\f. Finally the Hilbert-space Ti^ — on the two-torus has dimension N. Classical 
observables f{q,p) E C°°(T^) can be expanded in a Fourier series 

fiq,p)= Yl , (3.7) 

and their Weyl-quantisations are defined as 

Op(/):= 2: /.™T(^,^) . (3.8) 

Having reviewed the kinematic aspects of quantised torus maps, the next step is to quantise 
the dynamics, i.e. to find a suitable quantum time evolution operator U(i^a/3)- This will be 
achieved in terms of a Van Vleck-propagator Q with the generating function (p. 7]). However, in 
order to make this quantisation compatible with the discrete lattice, which the wave functions 
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are supported on, cf. (|3.3| ), one first has to find A'"-dependent rational approximations of a 
and /3, as has been pointed out by Marklof and Rudnick for quantum skew maps H. Following 
their approach, we approximate the irrational numbers by fractions with denominator A^, 
fulfilling 



a 



N 



< — — and 
- 2N 



N 



< 



1 

2iV 



(3.9) 



Thus, it is ensured that in the semiclassical limit — > oo we recover the irrational a and /3. 
With the generating function ( |2.7| ) and the approximations (|3.9| ) we now define the propagator 
in a mixed representation by 



{Q I ^{Kap) I 



1 



d'^Ga_±{q,p) 



dpdq 



exp 



9=^ 



2mNGa 6 

NN \ N N 



(3.10) 



with Q,P G Zat. Note that by ( |2.7| ) the pre- factor is constant. As can easily be verified 
using (|3.6|) we also have the identity 



u(ir, 



a/3 J 



A^ 



T O,- T -,0 



A^' 



(3.11) 



The time evolution operator U(/Cq,^) is a quantisation of the classical map ( |2.lD in the sense 
that in the semiclassical limit A^ — > oo it fulfils the Egorov property 



\V-\K^p) Opif)\j\K^p) - Op(/ o Kif,) 



(3.12) 



for all observables / G C°°(t2) and all fixed times t G IN. A proof of ( gl2|) IS given m 
appendix 

The propagator (|1i|) does not preserve all classical symmetries which were discussed in 
section ^. First of all the Hilbert space structure restricts the continuous families of translations 
in (|2.4|) and (^) to discrete translations 



A^' 



T — ,0 and T 0, — 



N 



(3.13) 



with c,d£ Ztv- But not all of these translation operators k := Op(7r) fulfil 

Tz\J{K^f3)K-^ = \J{K^f3) . (3.14) 

This is only valid for bc/i\i £ and '^'^Af G Z respectively. Consequently, only translations 



with D, 



b •= 



icd{b, N) 



h G 



and 



T O, 



ka 

a, 



with Da := gcd(a, A^) , ka G Zd, 



(3.15) 



(3.16) 
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are symmetries of the quantum propagator ( |3.11 ). Beside these quantum symmetries, which 



correspond directly to the classical translations ( |2.4| ) and ( |2.5| ), there are further quantum 
symmetries, 

f c d\ ^ he — ad ^ , 

which correspond to a combination of classical symmetries. Obviously, the last class of sym- 
metries ( p. 17 ) includes the previous ones. The inversion symmetry ( |2.3| ) is preserved and the 



corresponding operator I is given by 



m (Q) = i^i-Q) ■ (3.18) 

Summarising we emphasise that, although not all symmetries are preserved by the quantisa- 
tion of Kronecker maps, there are quantum symmetries in position and momentum, which are 
remnants of the classical translation invariance. However, it depends on the rational approx- 
imations for Q and /?, and thus also on the semiclassical parameter N, which of these phase 
space translations are conserved. Hence the characteristics of the quantum maps depend on 
number theoretical properties of a and /?. 



4. Trace formula and spectral statistics 

The purpose of this section is to investigate the spectrum of the propagator \J{Kai3)- The 
eigenangles i?j and the eigenstates \j) are defined by 



(4.1) 



Notice that by this definition the eigenangles are already unfolded, i. e. rescaled to mean density 
one. One could now explicitly calculate the eigenangles and the eigenvalues of the Kronecker 
map. However, we will not do so but instead discuss spectral properties in terms of exact 
and semiclassical trace formulae similar to the Gutzwiller trace formula [13|. This procedure 
will allow for an immediate generalisation to the case with spin which will be treated in the 
following sections. To this end recall that the spectral density can be expressed by traces of 
U'(/^„/3), 



(4.2) 



Thus, in order to derive a trace formula it is sufficient to calculate TrU^ One easily verifies 
by induction that the /-step propagator in mixed representation is given by, cf. ( 3.1C| ), 



al3) 



P 



exp 



b_ 

N N 



where 



G^lsiQ^P) =PQ + Wq - lap + 



a(3 



(4.3) 



(4.4) 
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generates the Z-step Kronecker map via 

and qo 



Pi 



dqi dpo 
In position representation the propagator reads 



Q 



al3) 



Q' 



exp 



((( - 1) 



ab 



},Q'+lamodN : 



and thus the trace is given by 
TrU'(i^, 



'('-!) ^fc 

TV e AT 2 OQ lamodNOo ibraodN 



(4.5) 



(4.6) 



(4.7) 



We only obtain a non-zero trace if both la and lb are multiples of A^. With gcd(a, N) = Da 
and gcd(6, A'') = Da we define Ma = ^/Oa and Mi, = ^/Dt respectively. Hence, in order to get 
get a non- vanishing trace, / has to be a multiple of both Ma and M^- 

The Kronecker map with translations "-/n and ^/n has a two-parameter family of periodic 
orbits 7. The periodic points are given by all {q,p) £ and the length of the periodic orbits 
is determined by 



= lcm(M„,M5) 



(4. 



Thus the condition for a non vanishing contribution in (^^) exactly singles out those values 
of / which correspond to multiples of l^. 

The exponent of ( [4. 71) has to be interpreted as the action of the periodic orbit in the 
following sense. In general one can pass from the generating function G^^(q',p) in mixed 
representation to a generating function 5'^^(g, g') in position representation by the Legendre 
transform 



pq 



Here, this step is only formal, since we cannot solve q' 

the initial momentum from (4.9). However, we can use ([4. 9]) to calculate the action 



(4.9) 

^'^"'dp'^'^^ ^'^^ P order to eliminate 



iJ a b 



Q + la Q 



N ' N 



b Q l{l-l) a b 



N N 



N N 



(4.10) 



of a periodic orbit, which appears in the trace formula for quantum maps. The last term 
involves the winding number lb in momentum and thus accounts for the difference of the 
action S on the torus and the action Si^^ , on the covering plane K,^. Therefore, 

the trace of is given by 



J V e do,l modMa 00,1 modMi, 



(4.11) 



Finally, we can express the spectral density ( |4.2| ) in terms of the periodic orbits of the classical 
map which we used to approximate the original map for quantisation. 



N N 



(4.12) 
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(a) Non-universal level spacing distribution of the 
Kronecker map. 



(b) Form factor of the Kronecker map which has 
peaks at r = fe h/N = ^3 with fe £ Z. 



Figure 1: Spectral statistics of a Kronecker map with a = (v^-i)/2, (3 = and N = 3000. 
The approximations ( |3.9| ) lead to a = 1854 and b = 2121. Thus we get Ma = 500 and 
Mb = 1000. 



Thus we have derived an exact trace formula for quantised Kronecker maps with fixed N. In 
the semiclassical limit the approximations a and b change and so does the length of the 
periodic orbit. For an irrational translation on the two-torus we have 

> 00 , (4.13) 

as should be expected for a uniquely ergodic map without periodic orbits. 

Prom the spectral density ( |4.12| ) it is clear that the Kronecker map does not show generic 
spectral statistics. Instead the behaviour for fixed N depends crucially on the approxima- 
tions (|3.9| ). This can be observed explicitly by investigating the form factor 

K{r) := 1 |Tr U^^ iK^(,)f - N6r,o = N6rN,o mod i, , (4.14) 



which has peaks if r is a multiple of h/N as shown in figure l(b)| . Also the probability density 
for the spacings sj := i^j+i — 'ffj, the level spacing density P{s), defined by 



behaves non-universal. In figure 1(a) the level spacing distribution 



I{s) := / P(s') ds' (4.16) 







is compared to the respective results of the circular orthogonal, unitary and symplectic ensem- 
bles (COE, CUE and CSE respectively) of RMT and the Poissonian level spacing distribution. 
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We only remark that in the Umit N ^ oo these statistics have no hmit distribution but depend 
crucially on number theoretical properties of the approximations (^^) , cf . the similar situation 
for the quantised skew map which was discussed in detail in Q. This behaviour is related to 
the (A^-dependent) quantum symmetries ( 3.17] ) which are remnants of the classical translation 
invariance (2^) and (|2.5| ). The connection becomes evident in the trace formula ( [4.12 ) where 
all degenerate periodic orbits collectively contribute to the spectral density. 



5. Kronecker maps with spin — generic statistics 

Having identified the quantum symmetries (and thus the classical symmetries) of the Kronecker 
map as being responsible for its non-generic spectral statistics, we now aim at breaking these 
symmetries without changing the characteristic properties of the classical system, i. e. without 
introducing any stronger chaotic property beyond (unique) ergodicity. To this end we couple 
the quantised map to a spin 1/2 |14,|l5|, i. e. we seek a quantum propagator acting on Tiiy 
instead of 7ij\f. The explicit coupling procedure can be motivated by a Pauli equation where 
the translational dynamics is periodically kicked in time [15|, yielding 



Uspin(i^a/3) = (U(i^„^) I2) Op (e--^('''^')) . (5.1) 



Here I2 is the 2x2 unit matrix, the components of B{q,p) are in C°°(T^) and cr denotes the 
vector of Pauli matrices with 

a. = (J J) > ^y= (i 0') = (i -1) ■ ^^-2) 

Thus, e^'^'^^'^'P^ takes values in SU(2) and is quantised by applying the Weyl quantisation ( p.^ ) 
to its components. The physical meaning of B{p, q) is that of a combination of magnetic fields 
and spin-orbit coupling terms, and thus \ispm{Kap) describes iterations of the quantum map 
alternating with spin precession. 

Note that the coupling of the spin degrees of freedom to the quantum map is in non- 
leading semiclassical order, as can be seen from (^.l|), where in contrast to ( 3.10| ) there is no 



pre- factor in the exponent of the spin part. Therefore the corresponding classical dynamics 
on the torus is not perturbed but only accompanied by a cocycle taking values in SU(2) as we 
will point out in the following. 

By multiple applications of the Egorov property ( |3.12| ) in the semiclassical limit we can 
write fl^ (by a slight abuse of notation we define \]{Kap) (g) I2 = \]{Kotp) ) 

Uspin(-f^a/3) 

\(tB 



\]{K^p) Op (e" 

U*(i^„^) U"(*-i)(i^„^) Op (e--^) U(*-i)(i^«^) U-(*-2)(K„^) . . . 

• • • \]\K^p) V-HK^p) Op (e---^) ViK^p) Op (e^'^) ^^"^^ 
V\K^f,) Op (e^-^-^ o K^-i) ... Op (e^-^-^ o K^f,) Op (e''^''^)) (l + O {N-') 



\j\K^fs) Op {dt) (1 + {N-^ 
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where the cocycle dt is defined by [|15|-|17| 

dt{q,p) := e'-^^^'/^^'P)) e^-^«/('?'f)) • • • e''^^^^'^') . (5.4) 
For later reference we remark that spin dynamics and classical map can be combined to form 



a skew product on x SU(2) by [|T8| 

Y : X SU(2) — > X SU(2) 

thus implying, Y^{q,p,g) = {K^p{q,p),dt{q,p)g). Notice that Y leaves the product measure 
/i = Hr^2 X hh invariant which consists of Lebesgue measure fir^2 on the torus and the normalised 
Haar measure fin on SU(2). Observables are now given by hermitian 2x2 matrices 



with f^^ G C°°(T^) which are quantised component-wise by (U). Prom (|1|) and Q now 
follows the Egorov property for Kronecker maps with spin, 

||U*pi,Op(F)U*pi,-OpK-i)Op(Foi^*^)Op(dt)||^^^^, =0(iV-i) . (5.7) 

for all fixed times t € Z. Clearly, when applying this Egorov property to a scalar observable, 
i. e. F{q,p) = f{q,p)t2, it reduces to ( ^.12 ). Thus the classical dynamics on the torus remains 



unchanged. 

We now present some numerical results on the spectral statistics of the Kronecker map with 
spin. For convenience instead of Op (^e^"'^(i'P)^ we choose Op (e"^"^^^''^) Op (e^*^"®^^^^) with 

(sin(27rg)\ /sin(27rp)\ 
sin(47rg) and B2{p) = sin(47rp) j . (5.8) 
sin(67rg) / \sin(67rp) / 

Due to the group property of SU(2) we can always find a unique B{q,p) with e"^"^^"^'^^ = 
gi(TBi((2) gi(TB2(p) g^j^jj ^]^g Moyal product of the Weyl quantisation, see, e. g., we have 

Op fe^'^-^(^'P)) = Op fe''^'^i('?)) Op fe''^-^2(P)) (l + O (N~^)) . (5.9) 



Therefore, the above semiclassical discussion applies to the numerical situation. Notice that 
the choice of the fields ( |5.8D breaks all symmetries of the quantum map \J{Ka/3). However, 
since the classical dynamics on the torus is not changed by the spin-coupling, being of higher 
order in A^, the quantum map VspiniKa/3) still corresponds to a non-mixing classical system. 
If the fields Bi and B2 generate a sufficiently random dynamics on SU(2) we may expect the 
skew product (|5.5| ) to be ergodic, but it can never be mixing since this would also imply the 
dynamics on the torus, given by to be mixing. 

In figure ^ we show the level spacing density and the level spacing distribution for a Kro- 
necker map with spin, observing good agreement with the CUE. This is consistent with the 



fact that Uspin(-f^o/3) has neither unitary nor anti-unitary symmetries. The form factor 2(b) 
shows typical fluctuations and when averaging (over an interval of length Ar = 0.2 with unit 
weight) one obtains again nice agreement with the CUE. 
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(a) Level spacing density in comparison to different 
circular ensembles. 



(b) In contrast to fig. 1(b) 
generic fluctuations. 



the form factor shows 



to 
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(c) The level spacing distribution I{s). 



(d) Smoothed form factor where we have averaged 
over At — 0.2. 



Figure 2: Spectral statistics of the Kronecker map with spin 1/2. As in fig. I we have chosen 
a = (V5-i)/2, p = V2/2 and N = 3000. Thus the Hilbert space has dimension 2N = 6000. We 
observe good agreement with the CUE. 



Thus we have provided an example in favour of the BGS conjecture, where only the weakest 
possible chaotic property, i. e. ergodicity, is fulfilled. The only difference to quantised Kronecker 
maps without spin, which show non-generic spectral statistics, is the absence of quantum 
mechanical degeneracies, which were discussed in the previous sections. In the following, we 
will explain that this is related to the lifting of degeneracies of periodic orbits in the semiclassical 
trace formula. We remark that for different values of a, /? G S^\Q we make similar observations. 
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6. Trace formula with spin 

The purpose of this section is to derive a trace formula for Kronecker maps with spin. Semi- 
classical trace formulae for flows with spin have been derived in [|^,^ and trace formulae for 
(perturbed) cat maps with spin have been investigated in |15|. 

First recall (|5.3D that the propagator of the Kronecker map with spin fulfils 

uU(^"/3) = U'(i^a/3)Op(dO (l + 0{N-')) , (6.1) 



where di denotes the cocycle (5^). Calculating the trace of (|6.lD thus reduces to deriving 
a (semiclassical) trace formula for U'(i^Q,^) with matrix elements of the quantum observable 
Op(/)) / = trc^i; where tr(-) denotes the trace on C^, i.e. we have 

Tr (U'(K,^) Op(/)) = I Q') {Q' I 0P(/) I Q) ■ (6-2) 

With the Weyl quantisation (p^) and the propagator (|4.6|) we obtain 



Tr(u'(K,^)Op(/)j =iV ^ /„™ew(— '('-i)«^)j„^,,^,diV<5,n,ifemodJV • (6.3) 

and thus get a non vanishing trace for 

n = la mod and m = lb mod N . (6-4) 

Let us discuss the first condition for fixed / 7^ in the semiclassical limit N ^ 00. We 
distinguish the following two cases: 

(i) la = mod A^: 

Then ( |6.4[ ) is clearly fulfilled for n = which leads to the same condition as without 
spin, see ( [1.11 ). The next values of n for which ( |6.4D holds are given by n = itA^. In 
consequence the smallest non-trivial value of n is given by min \n\ = N. 

(ii) la ^0 mod N: 

The smallest value of n for which ( |6.4[) is fulfilled is given by min |n| = min(Za, N — la). 
With the approximation (|3.9| ) we get min|n| = min(/A^a + 0(1), A^(l — /a) + 0(1)). Since 
I is fixed this leads to min|n| = 0{N). 

From (|) and (|i|) it is clear that in the semiclassical limit the first non-trivial n which yields a 
contribution to (|6.4[) is given by 



min|n| = 0{N) . (6.5) 



The second condition in ( |6.4| ) is discussed analogously. Hence in the semiclassical limit ( |6.3| ) 
reduces to the term with n = m = 0. All other terms are of higher semiclassical order since 
/ S C°°(T^) implies that the coefficients fnm of the Fourier series decrease exponentially. We 
therefore have 

TV (U'(i^a^) Op(/)) ~ iV/00 ew'(i-Oafe SoModNU So,imodM, , (6.6) 
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which can be written in terms of the periodic orbits of the approximated Kronecker map with 
length = lcm(Ma,M;,) and action from equation ( 4.10| ) yielding 



TV (uipiJK„^)) ~ iV ^ 6i,ki, e^-^S' / f ti di{q,p) dqdp , (6.7) 



T2 



where we have also substituted the observable f{q,p) by ivdi{q,p). As in the case without 
spin ( [4.11| ) we still have a two-parameter family of periodic orbits with action S'^'t. However, 
the contribution of each periodic point {q,p) G is weighted by the spin trace ii dki'y{q,p), 
thus breaking the degeneracies semiclassically without changing the periodic orbit structure of 
the classical torus map. 

Prom the Egorov property ( |5.7l ) it was already clear that the spin does not change the 
classical dynamics on the torus. This can now be nicely illustrated using the trace formula 
from above. To this end we evaluate the spectral density on a test-function gi;{L,'d), 



d,{L) := j gs{L,^)d{^) . (6.8) 

— oo 

Choosing 

(7.(L,i?) :=cos(t?L)e-'^'% (6.9) 



we obtain the so-called trace of the cosine-modulated heat-kernel |£0,21], and with the iden- 
tity ( |4.2| ) ds{L) becomes 

Therefore, by (|6.7]), in the semiclassical limit d^{L) has Gaussian peaks at multiples of the 
periods of periodic orbits. These can in turn be determined from the quantum eigenvalues. 
For a strongly chaotic map this is rather uninteresting since for each possible length L S IN 
an exponentially increasing number of orbits contributes to the amplitude. In our situation, 
however, there are only contributions for lengths L which are multiples of l^. Thus, each peak 
corresponds to a particular number k of repetitions of the family 7 of periodic orbits and we 
can observe how spin changes their effective multiplicity, i. e. the weight factor which these 
orbits contribute to the trace formula. 



In figure y we compare the cosine-modulated heat-kernel (|6.10D for the Kronecker map 
with and without spin 1/2. For both maps we have chosen the same parameters a and (3 as 
well as the same value of A^. Thus we get identical classical dynamics on the torus for the 
approximating maps which were used for quantisation. The dimension of the Hilbert space of 
the map with spin is twice the dimension of the Hilbert space of the map without spin. Thus 
the Heisenberg time T^™ for the map with spin is twice the Heisenberg time Th = for the 
map without spin and for the case with spin in ( |6.10| ) we simply have to replace A^ by 2N . 

For small L we get pronounced peaks at multiples of for both the map with and without 



spin. In agreement with the exact trace formula ( 4.12| ) for the pure Kronecker map without spin 



this remains true for all L. For ^sp\n{Kap), however, when increasing L we first observe that 
the height of the peaks decreases due to the spin weight ^2 ^"^di dqdp. Then, for L < T^'''", 
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(a) For lengths L <C T^^'" we get for both maps Gaussian peaks only for L = kl-y, G Z. The amplitude 
of the peaks differs for the map with and without spin 1/2 due to the contribution of the cocycle in the 
trace formula (6.7). 
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(b) If L becomes larger we get further peaks due to 
quantum fluctuations for the map with spin, but for 
the peaks at L = kl^ axe still dominating. 



(c) Finally, for lengths of the order rj^^'" one gets 
pure quantum fluctuations for the Kronecker map 
with spin. 



Figure 3: Cosine-modulated heat-kernel for a Kronecker map with and without spin. In both 
cases we have chosen = 3000. Thus the Heisenberg time is given by Th = 3000 and 
T^^^ = 6000 respectively. For both maps the approximations of a and /? yield Ma = 5 and 
Mf, = 2. Therefore the lengths of the periodic orbits are given by multiples of = 10. 

additional contributions appear which correspond to higher order corrections to the trace 
formula (|6.7|). Finally, for L ^ Tj^^™, these corrections dominate the behaviour. Summarising 
we have confirmed that in the semiclassical regime, i. e. for L <^ T^™, the coupling of the map 
to a spin only changes the amplitudes of the periodic orbit contributions but does not affect 
the classical dynamics on the torus. 

A similar discussion now applies to the spectral form factor, which for Uspin(-K^a/3) is defined 
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by 



K{t) :-- 



rpspm 



2N6r^ 



The smallest value of r for which K{t) does not vanish is given by 



I. 



lcm{Ma,Mh) 



rpspm 



^spm 



(6.11) 



(6.12) 



If a is irrational then Ma — > oo in the semiclassical limit — > oo [12, lemma 4.2], and similarly 
for Ml,. Thus if r tends to zero such that r < Tmin one concludes that 



lim K(t) 



N- 







»0 



(6.13) 



This is, however, the same result as we would have obtained in the case without spin, cf. sec- 



tion ^, which is consistent with both the non-generic behaviour in figure 1(b) and the generic 
behaviour in figure 2(d) . It is interesting to remark that with present methods, even in the 



regime which is usually dominated by the diagonal approximation |2^, 23|, it is not possible to 
distinguish semiclassically between generic and non-generic behaviour. This unusual observa- 
tion can be traced back to the fact that the Kronecker map has no periodic orbits and thus 
in the semiclassical limit, paradoxically, the periodic orbit contributions to the semiclassical 
trace formula do not determine the behaviour of K(t) in the vicinity of r = 0. Conversely, for 
small but non-zero r, in the semiclassical limit the form factor is dominated by corrections to 
the trace formula (J^^ which manifest themselves in figure |3| as additional peaks. 



7. Ergodic skew maps 

Another interesting class of parabolic maps on the two-torus are the so-called skew maps. 
We will consider skew maps of the form 

i2 , TpS 



5fl : — > T 



' (mod 1) . 



P 

qj yq + 2p 

For /5 G \ Q one has again a uniquely ergodic but non-mixing map, for details see [p^26|. 
If (3 is rational the map is in general almost integrable. This situation will be discussed in the 
following section, and is our main motivation for also analysing skew maps in addition to the 
simplest case of parabolic maps, namely translations of the torus. As in the case of Kronecker 
maps symmetries play an important role. Skew maps also have several symplectic symmetries. 
First of all there is a translation in momentum 



and a family of translations in position 

T,if)"Ll,) (-dl), 7^^S-. (7.3) 
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(a) The level spacing distribution for a skew 
map is shown. The distribution is obviously non 
generic. The level spacing density is given by 
P{s) = V3( S{s) + 5{s - 1) + 5(s - 2) ), see §. 



(b) The form factor has peaks for r = Yd — ^s, 
A; £ ISf. For multiples of the Heisenberg time Th = 
A'' one gets peaks of the height A*'. The average for 
r > 1 gives a mean degeneracy for the eigenangles 

of 5/3. 



Figure 4: Spectral statistics for a skew map with f3 = (v^-i)/2 and = 3000, implying D 
gcd(5, iV) = 3 and M = n/d = 1000. 



Moreover, skew maps are invariant under two anti-symplectic transformations, 



T : 



P 



-p + 6 



(mod 1) , 26 = 2/3 mod 1 , 



(7.4) 



Thus, in contrast to Kronecker maps, we now have a (non-conventional) time reversal symme- 
try. 

The first quantisation of skew maps, which fulfils an Egorov property for all G M, 
cf. (|3T^ ), was given by Marklof and Rudnik [^]. This quantisation is similar to the quantisa- 
tion of the Kronecker map which we presented in section |^. First of all one has to approximate 
the translation in momentum by a rational number f>/N with 



< 



2N 



b G 



(7.5) 



The propagator V{Sj3) can then be defined similarly to (|3.10| ) by using a generating function 
G{qn+i,Pn)- In mixed representation one obtains 



{Q\\J{Sp)\P) 



N \ N \n' N 



1 22n(^_p2_^Q(p+b)^ 



N 



e N 



and in momentum representation U(S'^) is simply given by 



[P\ViS0)\P') = e 



jp,P'+b ■ 



(7.6) 



(7.7) 



For this propagator the Egorov property was proven in [12|. Again the quantisation does not 
preserve all classical symmetries. For our purposes it is not necessary to give a complete list 
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(a) The level spacing distribution agrees with the 
CUE. 



1.2 ■ 




T 



(b) Smoothed form factor, where we averaged over 
At = 0.2. 



Figure 5: Spectral statistics for a skew map with spin 1/2. The parameters are identical to 
those chosen in figure |, i. e. N = 3000, /3 = (\^-i)/2, D = 3 and M = 1000. 



of all quantum symmetries, instead we only remark that there are always unitary quantum 
symmetries given by translations in position 

t(^-^,0^ , with D:=gcd{b,N) , keZo, (7.8) 

and an anti-unitary symmetry t := Op(r) involving a translation in momentum 



[T^](P)=v(-P+^j • (7.9) 

As for Kronecker maps one can calculate the eigenvalues and eigenvectors of the propa- 
gator ( |7.6| ) explicitly ||l^. Analysing the spectrum of skew maps one observes non-generic 
spectral statistics determined by the value of D defined in ( |7.8| ). Furthermore, there is even no 
limit distribution Q. The non-generic behaviour of the spectral statistics of a quantum skew 
map is illustrated in figure ^. For a semiclassical analysis we again use a trace formula. Details 
of the derivation can be found in appendix]^. With ( B.25 ), see also ( B.17| ), we obtain for the 
spectral density 



d(^) = ^ + lIL y y e^'^^'' e^-^^"^ ^ (7.10) 

^ ^ 27r 27r ^ ^ v/2ifcM 

with M := ^/d. The periods of periodic orbits of the approximating map are given by kM, 
and 1^ labels families of periodic points sharing the same action. These families form invariant 
manifolds lE^/afcA/ (V^)) see ( ^.21 ) below, which in particular are invariant under translations in 
position as is the classical map. Again we have high degeneracies in both the quantum spectrum 
and in the the trace formula which are the origin of the non-generic spectral statistics of the 
quantum map (see figure ^) . 
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As in the case of Kronecker maps we couple the quantised skew map to a spin 1/2 in order 
to break these symmetries. As we pointed out above the skew map also has symmetries in both 
position and momentum. Therefore, we use the same magnetic field B as for the Kronecker 
map ( |5.8| ), depending on both q and p. The same calculation as in section ^ then yields the 
Egorov property for skew maps with spin. 

In appendix |^ we also derive a semiclassical trace-formula for the skew map with spin. 
With ( B.25| ) the spectral density reads, 



(7.11) 



where now in contrast to ( 7.10 ) the different periodic orbits of one family, labelled by v, enter 
with different spin weights tr d. Hence the magnetic field again breaks the degeneracies in the 
periodic orbit sum ( |7.1l| ) , but does not change the classical dynamics on the torus in the sense 
that the set of periodic orbits remains the same and only their amplitudes are altered. 

In figure 5(a)| we show the level spacing distribution and the form factor for an ergodic 



skew map with spin 1/2. As in the case of the Kronecker map we observe good agreement with 
the CUE. Hence, we have a further example of a system which is only ergodic, but obeys the 
BGS conjecture. 

By using the same arguments as in section ^ and the semiclassical trace formula ( 7.11| ) we 
conclude that for the form factor of ergodic skew maps 

lim K{t) = (7.12) 

holds. But again, as in the case of Kronecker maps, the behaviour of K{t) for small but 
non-zero r cannot simply be determined by a semiclassical analysis. 



8. Almost integrable skew maps 

We now turn our attention to the interesting case of rational skew maps 



P + P\ , ... (8-1) 

with /3 E n Q, i. e. we can write [3 = B/M were M € IN and B E {0, 1, . . . , M - 1} are 
relatively prime. One immediately finds that invariant sets of Sp are given by M disjoint circles 

lEpo(/5):={(9,P)eT2 |p = po + |^, jGZm| . (8.2) 

If /? = then Epq(/3) = and thus the map is integrable. In all other cases one has M > 1 
and the map restricted to Epp(/3) can be represented as an interval-exchange transformation 
(for interval-exchange transformations see, e. g., [||,^,^ and references therein). Therefore, we 
will call this situation almost integrable or pseudointegrable. We only remark that for M — > 00 
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(a) Level spacing distribution for the map where 
one can see strong degeneracies in the spectrum. 



(b) Peaks in the form factor appear at k'^'/N = 
Veoo, with fc £ IN. The average for r > 1 gives a 
mean degeneracy of 13 for the eigenangles. 



Figure 6: Spectral statistics of an almost integrable skew map with N 
D = 600 and M = 5. 



3000 and (3 = s/s, i. e. 



we approach the situation of the last section, i.e. when increasing M, starting from M = 1, 
we observe a transition from integrable to ergodic. 

For the quantisation of 5/3 we can use the same procedure as in section |^ with the important 
difference that we will not approximate /? but instead restrict N to multiples of M, i. e. N = 
DM, D G M. As a consequence, the Egorov property, cf. (|3.12| ), is not merely an asymptotic 
relation but an identity. Following the notation of section |^ we will also write (3 = ^/n with 
b = BD. 

The spectral statistics of the skew map depends sensitively on the value of D, which ac- 
cording to (|7.8| ) determines the number of quantum symmetries. In Q it was shown that 
for /? ^ Q subsequences of quantisations with fixed D lead to different limit distributions as 
— > oo for various spectral statistics. In |12, corrolary 5.2] it was shown that the multiplicity 
of eigenphases of U(5/3) is bounded from above by c{e)D^/^^^ for each e > 0, where c(e) does 
not depend on D, i.e. for larger values of D one can obtain more degenerate spectra. Whereas 
for /3 E Q we have D = ^/m, i.e. D grows linearly with A^, for /? ^ Q it can be shown [Q] that 
as A ^ oo there exist subsequences with fixed D for any finite value of D. For example there 
is always a subsequence with D = 1. If /? ^ Q is badly approximable, D is bounded from above 
by c(e)A^/^^'^ for each e > lemma 4.2]. Thus, in general, for rational skew maps we obtain 
much larger values of D, which typically lead to higher degeneracies. This is nicely illustrated 
in figure 6(a)| , where /(O) is larger than in figure 4(a) . Roughly speaking, quantisations of skew 
maps with rational (3 are even more degenerate than those with /? irrational. 

We remark that one directly obtains an analytic expression for the non-generic form factor 
in figure S(b) from the trace formula for the skew map, cf. appendix ^ Moreover, the mean 
of the form factor for r S> 1 gives the average of the degeneracies, which now is 13 in contrast 
to ^or the ergodic example in the previous section. 

When now coupling a spin 1/2 to the map as in the previous cases, we again expect that 
this procedure lifts the degeneracies of periodic orbits in the trace formula and the degenera- 
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(a) The level spacing distribution shows intermedi- 
ate statistics. 



(b) The smoothed form factor also shows interme- 
diate statistics. For t = one observes a value 
between and 1. 



Figure 7: Spectral statistics of an almost integrable skew map with spin 1/2. We have chosen 
the same parameters as in figure ^. Therefore we have A'^ = 3000, (3 = 3/5, D = 600 and M = 5. 



cies of quantum eigenphases and thus restores generic behaviour. The resulting level spacing 
distribution /(s) and the averaged form factor K{t) are shown in figure |^. The non-generic 
features in figure ^ caused by the high degeneracies, have vanished and we obtain an level 
spacing distribution between COE and Poisson showing level repulsion on short scales but for 
larger s staying below the RMT curves. The form factor starts with a non-zero value at r = 
and then approaches 1. 

The intermediate level spacing distribution 
of almost integrable - or pseudointegrable - 
systems, cf., e.g., |2^, 2^, is often associated 
with so-called Semi- Poisson statistics [30|, 31| 
which yields 



Isp{s) 



4s e" 



-2s 



and 



1 - (2s + l)e 



-2s 



i.3) 



At a first glance this appears to be consistent 
with figure |^. However, when looking more 
carefully, in figure ^ we still observe fluctua- 
tions about /sp('5) when only slightly chang- 
ing N. The value of K{0) also changes and, 
thus, from our data we cannot judge whether 
the statistics will converge to a limit distribu- 
tion as — > 00. However, considering our 



Semi-Poisson 
Poisson 
COE 




Figure 8: Level spacing density of the almost 
integrable skew map with spin in comparison 
to Psp(s). 



previous discussions it is not clear whether one should expect convergence. We have argued 
that in order to observe generic statistics one has to break the degeneracies of the quantum 
skew map. Since these degeneracies grow with D this can be nicely achieved for irrational skew 
maps where we obtain CUE statistics. For rational skew maps, on the other hand, the quickly 
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(a) The level spacing distribution varies with A'^. 
We get fluctuating about the Semi-Poissonian dis- 
tribution. 



(b) The form factor also varies and even the value 
for T = changes slightly with A'^. 



Figure 9: Comparison of spectral statistics for an almost integrable skew map at different values 
of A^. We show the level spacing distribution and the smoothed form factor where we averaged 
over At = 0.2. 



growing D might prevent us from restoring genericity and we might only observe a tendency 
towards Semi-Poissonian statistics. 



However, based on the semiclassical trace formula ( B.25 ) for the skew map with spin we 
can give a semi-heuristic argument yielding the value to which K{{)) should converge in the 



semiclassical limit. If we insert ( B.25 ) in the definition of the form factor 



Kir) = — 
^ ' 2N 



2N6oi , T = ^ 



I 

2N 



(8.4) 



and perform the diagonal approximation |22,23], i. e. in the double sum over families of periodic 
orbits, labelled by i^, we keep only terms with like actions, we obtain 



'tvdi{q,p)) 



(8.5) 



where 



- y 



M 



(8.6) 



denotes the average over the invariant manifold Ep{^)(/3). 



Due to the anti-symplectic sym- 
metries (7.4) there are additional degeneracies in the actions. However, assuming that the 
corresponding spin weights trd; are sufficiently uncorrelated we do not take correlations be- 
tween these terms into account. Thus there appears no additional factor of 4 in ( |8.5[ ). In order 
to determine K[Q) we average K{t) over an interval (0, r*) and simultaneously perform the 
limit N ^ oo and r* ^ such that t*N oo. Since then the dominating contributions come 
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from large I we are interested in the behaviour of ( ^.51) as / — > oo. In this hmit we can replace 
the sum over v by an integral obtaining 

(8-7) 

fcez 

with the combined average 

1/a/ 

ai:=M j [(trdKg,p)>]E^,(^)]' dp'. (8.8) 



Performing the r-average or /-average, respectively, as described above, yields 
1 1 ^ 

E:(0)~- hm - Y^akM. (8.9) 

k=l 

and, since due to | tr d/| < 2 we have < < 4, the value of K(0) is restricted to the interval 
[0,2]. 

We will now give a heuristic argument why it is reasonable to expect K{0) ~ 1/2. To this 
end consider [{tr di{q,p))^ S™^^ skew map Sf^ restricted to Ep/(/3) is ergodic for 

almost all p' we can replace ( • ,(^) by some time average ~* along an orbit with initial 
condition {qo,Po) £ Ep/(/?). We now assume that for large / the spin weights ti di{Sp{qo,po)) 
are uncorrelated, which for the magnetic fields ( |5.8D is a reasonable assumption. Therefore, we 
have (trd; )^ = (trd^)^ , which again due to the ergodicity of Sfs on Ep'(/3), implies 

ai=M [ {{iTdi{q,p)f)^ dp' . (8.10) 







If we also assume that the skew product : {q,p,g) 1— > {S^i3{q,p),di{q,p)g), cf. (|5.5| ), is 
ergodic on Ep'(/?) x SU(2) for almost all p' (which, again referring to the fields {f>.8\), is a good 
assumption) we obtain, cf. jlS, 18|, 



1 ^ 1 ^ f 

K— >oo A ■'^^ — ' L— ►cxD ■^^ — ' J 
k=l 1=1 

where denotes the normalised Haar measure on SU(2). Thus we find K{{S) ~ 1/2, which is 
in agreement with the Semi-Poisson distribution |]30| , |3ll| . However one should note that both 
in the light of our previous discussion about the slow numerical convergence and considering 
the various assumptions made in the last paragraph this should not be viewed as a derivation 
but rather as a semiclassical illustration of our numerical findings. A final decision about 
whether the quantised skew map with spin shows Semi-Poissonian statistics requires further 
investigations, both numerical and semiclassical. 
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9. Conclusions 



We have investigated quantised parabolic maps on the torus which have highly degenerate 
spectra showing non-generic spectral statistics. An outstanding property of their classical 
limit is the high degeneracy of periodic orbits occuring in continuous families. For ergodic 
maps we have shown numerically that both the classical and quantum degeneracies can be 
lifted by coupling the map to spin 1/2, thus leading to generic spectral statistics following RMT 
predictions. Our numerical findings are supported by a semiclassical analysis. Hence we have 
given evidence that in generic situations ergodicity might be a sufficiently strong condition in 
order to observe RMT statistics, sheding new light on the weakest possible requirements for 
the BGS conjecture. 

Moreover, we have investigated an almost integrable map which also has a highly degenerate 
spectrum. Again by coupling the map to a spin 1/2 we numerically observed a transition to in- 
termediate statistics. We have also given a preliminary semiclassical discussion on whether this 
situation can be described by the so-called Semi-Poisson statistics. However, a final decision 
on this particular issue is beyond the scope of the present paper. 
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A. Egorov theorem for the Kronecker map 



For the sake of completeness we give a proof of the Egorov property ( p. 12 ) for the quantised 
Kronecker map, cf. |12|. To this end first consider the Weyl quantisation of a time evolved 
observable f{q,p), 



Op(/oi^*^) = Op Yl 



(m((?+o)-n(p+/3)) 



(A.l) 



Ef 2nit{ma-n/3) rp f ]]_ 



with t € Z. On the other hand, by virtue of the commutation relations ( |3.lD , for the quantum 
time evolution we have 

U-(i^.,)Op(/)U*(i^.,)= Y /.™e-*^T(^,^) . 



(A.2) 



Thus, the discrepancy between classical and quantum mechanical time evolution is given by 



A fit) := \\V-'iK^p) Op{f)V'{K^p) - Op(/ o 



n m 
iV' N 



(A.3) 



< 



0^;+ ma — r 

g27nt^^ 



^2TTit{ma—nf3) 



^ n m 



23 



With ||T(;^ 



A" N) WUt 



1, and an application of 



jgizx _ gizyj ^ l^ll^ _ ^1 ^ X^y^Z^'R, 

we have 



(A.4) 



/\f{t)<2At\ (|- 



a 



n + 



n\ \fn 



(M) 



-^7r|t| ^ (|n| + |m|)|/„ 



(A.5) 



which concludes the proof, since / G C°°(T ) implies that fnm decays exponentially as n, m 
oo, 



B. Trace formula for skew maps 

Before calculating the trace formula for the skew map we will summarise the structure of 
periodic orbits 7 for a skew map with a rational translation ^/n. From the definition of the 
map we get the condition for periodic orbits 



qo + 2lpo+Kl-l)^ 



' 'Po 

.10 



(mod 1) 



The period is determined by the condition in momentum. With 



D := gcd(6, N) and M :-- 



N 

'd 



we get 



= kM , with k £M . 
The other equation gives the momenta of the periodic orbits, 

V 



Po 



2kM 



with u E 7,2kM 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



There is no restriction on the initial position go and thus the periodic orbits occur in one- 
parameter families. 

After this preliminary remark we now turn to the calculation of the trace formula. Starting 
point is the one-step propagator ( f7.6| ) which in position representation reads 



Q 



u(^/3)|Q') = i E 



exp 



(B.5) 



One easily verifies by induction that the /-step propagator is given by 
{Q\v\Sf,)\Q') = ^ exp 



(B.6) 
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where 



Cl/^ {qi,po) ■= qi (po + ) - ^(po + ) 

^ ^ A=0 ^ ^ 



(B.7) 



denotes the generating function for the /-step skew map from (pQ, qq) to (pi, qi), i. e. 

Pi = K and go = ^ • (Bi 

dqi dpo 

Applying the Poisson summation formula to the propagator ( |B.6| ) yields 

1 



Q 



\j'iSf3)\Q') = Y, /exp 



Q \_q;_ 



27TiN[Gi,^ { -^p + np 



dp . 



(B.9) 



We can split the sum over n by writing n = 2lfi + i^, fj, 1i, u £ and thus we get 



27riiV ( G^i, l9.^p\-9^p+ (21 Pi + v)p 



With the substitution p' = p + n we can rearrange the integrals such that 

00 

(Q I U'(5,) \Q')=Y: I exp [2vriiV (gI/, " (| " P' 

The (Gaussian) integral can be calculated yielding 

u'(5/3) Q' 



dp . 
(B.IO) 



dp' . (B.ll) 



Q 



21N 



h E ^^p 



ueZ2t 



2TTiNS: 



■/^ \ AT + AT 



(B.12) 



where we have taken the freedom to change u to —u. We denote by Sl^j^{q,q') the Legendre 
transform of G[, {quPo) with respect to pq, 



b/N 



2 l+l b l-l b 



(B.13) 



A=0 



with 



Pi 



dSi/^ {qi,qo) 
dqi 



and pq 



dSi/^ iqi,qo) 



dqo 



(B.14) 



Therefore, the initial momentum p^'^^ of an orbit from Q'/n to Q/n + v with length I is given 
by 



p 



Q-Q' V l-l b 

+ 7TT + 



21N 21 2 N 



(mod 1) 



(B.15) 
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where v labels the 21 different orbits from Q' jN to QjN. Hence the propagator in position 
representation is given by the sum over orbits from Q'/iV to QjN . We can now calculate the 
trace of the propagator (|B.ll ) , 



TrU'(5^)= ^ (Q \J\S(,) 



Q 



(B.16) 



By using eq. ( B.13| ) we get 

77 



TrU'(S^) 



^^«,feM ^ exp 



27riiV§ 



with 



Sb/N {q + z^, 



(B.17) 



(B.18) 



One easily verifies that the right-hand side of ( [B.18| ) does not depend on the position q. Notice 
that I^/n is the winding number in momentum. Thus for periodic orbits with fixed momentum 
p we have the same action S[, for all g G S^. 

It remains to show that the structure of the trace formula remains the same when the 
coupling to the spin degrees of freedom is introduced. To this end we have to calculate (cf. 
section ^ and |jl5| ) 



1VUipi„(5^) = Tr(U(S^)0p(/)) 
with f = trdi. With (U) and j ^J^ ) we have 

ji 



(B.19) 



. nm 

X fnm exp ( ITT— + 27ri ^ 



N 



(B.20) 



For / G C°°(T^) we may expand the action as N ^ oo (cf. p5[| ) 



S\ 



Q 

b/N \ ^ 



Q — m 
N 



s; 



Q Q 



b/N \jy 



N 



dqo \N 



u, ■ 



Q\ m 



N N 



§1 +l^9. + !IL„M 



(B.21) 



Now by the same considerations as in section ^ for / fixed and N large the sum over the 
positions Q becomes 



E exp (^(/6 + n)Qj ~5„o^iV5i,fc 



kM 



(B.22) 



i. e. it selects the periods kM of periodic orbits of the skew map. Thus, the trace of the 
propagator reads 



IM 




(B.23) 
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By ( |3.7D the sum over m is given by f{p^'^\q) dq, i.e. resubstituting f = tTdi we obtain 

1 

Observing that = if p^"'^ G Ep(.)(Vw), cf. which can be checked by direct com- 

putation, instead of summing over the momenta p^'^^ of periodic points we can sum over the 
invariant manifolds Wipg(}/N) and their disconnected components. Finally we have derived the 
trace formula 

1 

where, p^^^ S Ep{i,)(yAf) can be chosen such that < p^'^^ < i/a/. Now the right-hand side 
of ( |B.25| ) has to be interpreted as follows: We have obtained a sum over all periodic orbits of 
the (approximating) classical skew map. The periodic orbits appear in families labelled by u 
and the spin weights trd/ are integrated over the corresponding invariant manifold Ep{^)(*/Af). 
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